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Abstract 

We study quantum phases of ultracold bosonic atoms in a two-dimensional optical su- 
perlattice. The extended Bose-Hubbard model derived from the system of ultracold bosonic 
atoms in an optical superlattice is solved numerically with Gutzwiller approach. We find that 
the modulated superfluid(MS), Mott-insulator (MI) and density-wave(DW) phases appear 
in some regimes of parameters. The experimental detection of the first order correlations 
and the second order correlations of different quantum phases with time-of-flight and noise- 
correlation techniques is proposed. 
PACS: 03.75.Lm, 03.75.Nt, 32.80.Pj 
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The system of ultracold atoms in optical lattices, which provides an intriguing environment to 
study strongly correlated condensed matter systems and quantum information due to the com- 
plete control over the system parameters, has been extensively investigated theoretically and ex- 
perimentally. The superfluid-Mott-insulator transition of ultracold bosons in optical lattices has 
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been theoretically analyzed and experimentally demonstrated 

tum liquids with atomic gases has achieved 3, 4|. A Bose glass phase has been observed in an 
optical latticejSi]- Repulsively bound states in an optical lattice have been observed in a recent 
experiment Q]. Atomic Bose- Fermi mixtures in an optical lattice have been studied 7j_8|. Ultracold 
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fermionic atoms in optical lattices have also been investigated recently [9|, 
quantum phase transition, nonlinear dynamics of a dipolar Bose-Einstein condensate, phase dia- 
gram of two-species Bose-Einstein condensates and magnetic soliton and soliton collisions of spinor 
Bose-Einstein condensates in an optical lattice have been studied 
of quantum information processing in optical lattices have been proposedj. 
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The quantum 

phases of ultracold bosons in one-dimensional optical superlattices have been discussed in some 
literatures 
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261 ] ■ However, the system of ultracold bosons in a two-dimensional 



optical superlattice has not investigated as yet. 



* Electronic address; jmhou@seu.edu.cn 



2 



In this paper, we will study the system of ultracold bosonic atoms in a two-dimensional optical 
superlattice. We will show that novel quantum phases, such as modulated superfluid, density wave, 
and Mott-insulator phases, appear in some regimes of quantum parameters. We will display the 
first order correlation function and the second order correlation function in the lattice momentum 
space for different quantum phases in order to identify them in experiments with the time-of-flight 
and noise-correlation techniques. 

Ultracold bosons with repulsive interactions in an external potential can be described by a 
second quantized Hamiltonian as follows, 

H= [ d^ripUr) - — V^ + V(r)-fl -^(r) + - f d^rijUrUUrU(rU(r) (1) 
J 2m 2 J 

where '(/'(r) and iP'^{y) are the boson field operators for annihilating and creating an atom at site r 
respectively; g = Airasfi^ /m, where Ug is the s-wave scattering length, which is positive for atoms 
with repulsive interactions; m is the mass of atoms and p, is the chemical potential; V{r) is an 
external potential. In our model, the external potential V{r) is a two-dimensional superlattice po- 
tential such as V{x,y) = Vi cos^(/cx) cos^(A;y) -|- V2[cos^(A;x/2) cos^(/cy/2) -|- sin^(A;x/2) sin^(A;y/2)], 
which is schematically displayed in FIGlH This superlattice potential consists of two sublattices 
that we denote as sublattice 'A' and sublattice 'B' with the potential depths Va and Vb, respec- 
tively. Such a two-dimensional superlattice can be built with optical techniques. For simplicity, we 
do not consider the shallow harmonic trapping potential, which is usually required in experiments 
to trap atomic gases. 

In periodic potentials, the atomic states are represented by Bloch wave functions, which can 
be expanded in terms of a set of Wannier states that are well localized at lattice sites. Because of 
the translational symmetry of each sublattice, the Wannier function Wi{r) is equal to WAi^ — R-O 
for sublattice A and is equal to tf^lr — Rj) for sublattice B, where Rj is a lattice vector. At zero 
temperature, the atoms are in the lowest band states. Expanding the field operators in the lowest 
vibrational Wannier states as ijj{r) = J2ibiWi{r), Eq.([T]) reduces to the extended Bose-Hubbard 
Hamiltonian, 

H = -*E[^k,i,)W,%) + ^L,.,)%.,%+i) + H-c.] -E^«^* + ^Ef^^"*(^^-i) (2) 

i i i 

where {ix^iy) is the coordinate position of lattice site i in the two-dimensional optical superlattice. 
^iix,iy) ^^'^ i ) represent the operators for annihilating and creating an atom at lattice site 
i = (ixjiy) respectively; t = J dr w\{r — rj)[— -|- y(r)]iz;B(r — Tj) is the tunnelling parameter 
between adjacent lattice sites in the different sublattices A and B. Ui is the repulsive on-site 
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interaction between atoms at lattice site i; /ij is the effective chemical potential corresponding 
to lattice site i. In the two different sublattices, the effective chemical potential and the on-site 
interaction are 



where ^ = p. — J drV{r)\wA{''c)\'^ and S = J drV{r){\wB{Tc)\'^ — |iL'yi(r)p); Ua = 9 J dr\wA{T^)\^ /rn 
and Ub = g J dr\wB{Y)\'^ /m are the on-site interactions on the two different kinds of sublattice 
sites. 

With Gutzwiller's ansatz, we can write the wavefunctions of the system as[l|, 



where \n)i is the Fock state with n particles at the lattice site i and /A are the variational 
parameters, which satisfy the normalized condition J2n = 1 foi' each lattice site i. To find the 

ground state of the system, we minimize the expectation value of the Hamiltonianfp]) . {^g\H\^g): 
with the variational parameters {fn^}, under a given chemical potential /x. In our work, we use 
conjugate gradient algorithm to minimize the expectation value of the Hamiltonian. The system 
contains 10 x 10 lattice sites and the truncated number of Fock states is 6. To eliminate the 
boundary effects, periodic boundary conditions are used. 

To identify different quantum phases, we define the static structure factor as follows: 



where Rj is a lattice vector and N is the number of lattice sites. The diagonal long-range order and 
the off-diagonal long-range order in the system of ultracold bosonic atoms in a two-dimensional 
optical superlattice are measured by S{-k,tt) and (j) = (6), respectively. The Mott-insulator phase is 
a state with particle number density pinned at an integer such as n = 1, 2, • • •, which corresponds to 
a commensurate filling of the lattice. The Mott-insulator phase is characterized by ^(Tr, vr) = and 
<j) = 0. In the density wave phase, the particle number density is modulated by with the period 
of twice lattice constant and is pinned at two different integers in the two different sublattices. 
The density wave phase is a superlattice crystal with the diagonal long-range order measured by 
S{tt,tt), and without the off-diagonal long-range order. Therefore, ^(Tr,?!) / and cp = feature 
the density wave phase. In the modulated superfluid phase, the diagonal long-rang order and 
the off-diagonal long-range order coexist, so the modulated superfluid phase is characterized by 
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^(Tr, vr) 7^ and </> 7^ 0. The conventional superfluid phase, which is absent in our model, is 
characterized by ^(vr, vr) = and 7^ 0. 

FIG. [2] shows the phase diagram of ultracold bosonic atoms in an optical super lattice for Ub = 
0.5Ua and 6 = 0.2Ua, which is obtained by solving the Hamiltonian ([2|) numerically with Gutzwiller 
approach. This phase diagram is a result of competition among the tunnelling term and the two 
different on-site interactions in the two different sublattices. From this figure, we can find that Mott- 
insulator, density wave and modulated superfluid phases appear in some regimes of parameters. 
The Mott-insulator and density wave phases are denoted as {nA,nB} for their filling numbers in 
the sublattice A and B. The second and fourth lobes counted from the bottom, denoted as {1, 1} 
and {2,2} respectively, belong to the Mott-insulator phase. For the Mott-insulator phase {1,1}, 
the critical value of t/UA for the MI-MS phase transition is 0.031, while the corresponding value 
of the MI-SF transition in a two-dimensional optical lattice is 0.043[l|]. Thus, the Mott-insulator 
phase is more instable in optical superlattice than in optical lattice due to competition between 
the different on-site interactions in the different sublattice in the optical superlattice. The first, 
third, fifth and sixth lobes, denoted as {1,0}, {1,2}, {2,3} and {2,4} respectively, belong to the 
density wave phase. The reason that the density wave phase {1,0} appears is the difference of 
the effective chemical potentials for the two different sublattice, while the other density wave lobes 
are the result of competition among the tunnelling term and the two on-site interaction in the 
two sublattices. The other regime except the Mott-insulator and density wave lobes, where the 
diagonal long-range order and the off-diagonal long-range order appear simultaneously, belongs to 
the modulated superfluid phase. In our model, the superfluid phase, which conventionally appears 
in the optical lattice-cold atom system[l|, Q], is modulated into the modulated superfluid phase by 
the periodic modulation of the optical superlattice. 

The techniques to detect the above quantum phases are accessible for now. Altman et al. 
have proposed to utilize shot-noise correlations to probe complex many-body states of trapped 
ultracold atoms, e.g., Mott statesi!27]. Scarola et al. have suggested to use the noise-correlation 
technique to detect the supersolid phase of cold atoms with long-range dipolar interactions in an 
optical lattice[28]. The correlations for the case of a bosonic Mott insulating state jisl], and the 
pair-correlated atoms created by dissociating weakly bound diatomic molecules near a Feshbach 
resonance [29|] have recently been observed through shot-noise correlations techniques experimen- 
tally. 

For the time of flight detection, the atoms interact weakly after the trapped optical lattice is 
adiabatically turned off, and the number of atoms at position r in the expanding cloud after a 
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time T is given by {n{r))T ~ |t()(k(r))p/9o(k(r)), where w is the Fourier transform of the Wannier 



function 



27 



28| . The relation between the lattice momentum and the position in the expanding 
cloud after an expanding time T is k(r) = aomr/hT, where m is the mass of atom and oq is the 
width of the Wannier state in the lattice, po is a first order correlation function in the lattice 
momentum space which is defined as follows, 

poik)=J2e'''-^'''-''^Hblb,), (5) 

id 

where Rj is a lattice vector. The first order normalized correlation function is defined as 
p(k) = pq(]<.)/N, where is the number of lattice sites. The shot-noise correlation function in the 
expanding atomic cloud after an expanding time T is ^(r, r') = (n(r)n(r'))r — (n(r))2^(n(r'))7^, 
which is proportional to momentum correlations in the ground state of the optical lattice trapped 
system, 

e(r,r') ~ |tI;(k(r))|2|u;(k(r'))pGo(k(r),k(r')), (6) 

where Go(k, k') is the second order correlation function in the lattice momentum space, the second 
order normalized correlation function is defined by G(k, k') = Go(k, k')//9o(k)/9o(k'), which can be 



written as 



28| 



^ '^ Po(k)po(kO hT Iwik'Wpoik'Y 

where the 6(k — k') term is the autocorrelation term and will be dropped in the subsequent 
discussion because it only contributes to the second order correlation function for k = k' and 
does not contain the characterizing information about the correlation function distribution in the 
momentum space. 

To identify the quantum phases, such as the Mott-insulator, density wave and modulated super- 
fluid phases, in experiments, we pick out three points a, (3 and 7, representing the three quantum 
phases respectively, from the phase diagram shown in FIGEJ For the three picked points, we cal- 
culate the distributions of p(k) and G(k, 0) in the lattice momentum space, which are shown in 
FIG. El The first order normalized correlation function is almost a constant for the Mott-insulator 
and density wave phases, while, for the modulated superfluid phase, it has peaks at points (0,0), 
(0,ib27r), (ib27r, 0) and (ib27r, ib27r) in the lattice momentum space, which indicate that the off- 
diagonal long-range order exists. The second order normalized correlation function has peaks at 
points (0,0), (0,ib27r), (ib27r,0) and (ib27r, ib27r) in the lattice momentum space for the Mott- 
insulator, while, for the density wave and modulated superfluid phases, it has additional peaks at 
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(±7r, diTr), which indicate that the diagonal long-range order exists. These characterizing features 
of correlation functions in the lattice momentum space can be detected by probing the correlations 
in the coordinate space of the atomic cloud after an expanding time T with the time-of-flight and 
noise-correlation techniques. 

In our work, for simplicity, wc do not consider the shallow harmonic trapping potential, which 
is usually required to trap atomic gases in the practical experiments. In fact, the presence of a 
shallow harmonic trapping potential does not change the intrinsic features that we have discussed. 
The energy offset due to the shallow harmonic trapping potential can be combined in the effective 
chemical potential. Thus, the effective chemical potential varies with the position, so the phase 
separation may occurs in the practical experiments. For the time-of-flight and noise-correlation 
detections, in the practical experiments where the shallow harmonic trapping potential exists, we 
may obtain lower peaks at the edge points than at the center in the expanding atomic cloud. 

In summary, we have investigated the quantum phases of ultracold bosonic atoms in a two- 
dimensional optical supcrlattice. First, we derived the extended Bose-Hubbard Hamiltonian with- 
out nearest-neighbor interactions. In our model, through solving the extended Bose-Hubbard 
Hamiltonian with Gutzwiller approach we found that the modulated superfiuid, density wave, and 
Mott-insulator phases appear. The detection of these quantum phases with the time-of-fiight and 
noise-correlation techniques was proposed. We have calculated the first order and second order 
correlation functions for different quantum phases and discussed their differences to distinguish 
them in experiments. 
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FIG. 1: (color online), (a) Schematic diagram of a two-dimensional optical superlattice potential, where x 
and y are coordinates in units of lattice costant a and V{x^y) represents the potential energy at position 
{x,y). (b) The profile of the optical superlattice potential shown in (a) at y = 0. Here Va and Vb represent 
the lattice depths of the sublattices A and B respectively. 
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FIG. 2: Phase diagram of ultracold bosonic atoms in an optical supeiiattice for Ub — 0-5Ua and 5 — 0.2Ua- 
t/UA and ij./Ua are the dimensionless hopping matrix element and the chemical potential respectively. In 
this diagram, MI, DW and MS denote the Mott-insulator, density wave and modulated superfluid phases 
respectively, {ua, ns} represents the average particle density in sublattices A and B when the system is in 
MI or DW phase. The points marked with 'x', of which the first order correlation function and the second 
order correlation function are shown in FIGlHl are denoted by 'a', '/?', and '7'. 
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FIG. 3: (color online). The first order and the second order normalized density correlations of different 
quantum phases as functions in the lattice momentum space, (a) and (b) are respectively the first order 
and the second order normalized density correlations in the lattice momentum space of a Mott-insulator 
phase at the point denoted by 'a' in FIGlH (c) and (d) are respectively the first order and the second 
order normalized density correlations in the lattice momentum space of a density wave phase at the point 
denoted by '/?' in FIGO (e) and (f) are respectively the first order and the second order normalized density 
correlations in the lattice momentum space of a modulated superfluid phase at the point denoted by '7' in 
FIGH 



